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I. THE BASIC IDEA 
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The "unparticle" proposal [l[ is based on the assumption that there is a type of new physics (NP) with the peculiar 
properties of being asymptotically free in the ultraviolet (UV) and conformally invariant in the infrared (IR). It is 
assumed that the NP interacts weakly with the standard model (SM) through the exchange of a set of heavy mediators 
of mass Mu . Though not explicitly stated it is also tacitly assumed that the theory will have an ultraviolet completion, 
whose details (aside form the existence of the above-mentioned mediators) are left unspecified. The NP sector has two 
relevant high-energy scales: Mu and Ay the scale at which conformal invariance sets in. One must have Mu > Au and 
we will also assume Au > v = 246GeV, the electroweak scale. The basic example for this type of NP is provided by 
a model proposed by Banks and Zaks Qj. In the following we will denote the UV phase of the NP as the Banks-Zaks 
(BZ) phase, and the IR phase as the unparticle (U) phase. 

At energies below Mu all mediator effects are virtual and generate effective interactions between the new physics 
and the Standard Model (SM) sectors. The dominant interactions are assumed to be of the form 



C(UV) = cbzM^OsmObz 



(1) 



where Osm is a local, gauge-invariant operator constructed out of the SM fields and their derivatives and, similarly, 
Obz denotes a local operator constructed of the NP fields in the BZ phase. In the infrared these type of effective 
interactions suffer from strong renormalization effects leading to the replacement Oqz — > A^f z ~ du Ou, where dgz 
and du denote the scaling dimensions of the operators Ogz and Ou respectively. We then have 



C(IR) = cuM-"0 
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(2) 



The construction of the operator Ou is in general a difficult task, and requires detailed knowledge of the NP theory. 
However, for most calculations this is not needed since one is interested in operator correlators and these are very 
strongly constrained by conformal invariance. If one considers the effects of a single operator Ou then for the purposes 
of calculating cross sections and for the effects on standard cosmology one only needs the following density of states, 



2\d u 



d i xe^([Ou(x),Ou(0)]) = A du 6{q'){q') 
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(3) 



where A n = (47r) 3 ~ 2 "/[2r(n)r(n — 1)] and (• • •) denotes the thermal average at temperature 1//3, and where we 
assumed Ou has bosonic character. For collider applications one simply lets (3 — > oo. Using this expression one can 
determine the effects of this type of NP on various experimentally interesting processes in terms of a few parameters 
(k, du, Mu and A^). In the following we will study some of the effects of this type of NP on standard cosmology 
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II. THERMODYNAMICS 



The conformal invariance requirement amounts to the assumption that the NP beta function vanishes at g = g* 7^ 0. 
In this case the (3 function and g will behave qualitatively as in Fig. [1] 
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FIG. 1: Renormalization group behavior of unparticle theories. 

For such theories the trace of the energy-momentum tensor Q obeys 

K) =Pu- iPu = W2d) (N [FF F a fiv\) = (P/2g)bT i+s 



(4) 



(N denotes normal ordering). It follows that (9£) will vanish in the IR (since it is oc (3) and we expect pu — 5Pu at 
low temperatures. The leading corrections are produced by (N [F£"F a as indicated above, with 6 the anomalous 
dimension of this operator. 

Using then standard thermodynamic relations we find 



p u = <jT 4 + A(l + 3/S)T 



4+5 . 



Pu = <jT 4 /3 + (A/8)T 



4+5 



(5) 



valid in the IR regime. In the UV limit, we will also have /?np oc T 4 (up to logarithmic corrections) since the theory 
is asymptotically free. Then 



P = ~~2 



5NP 



gsz phase 
gu U phase 



(G) 



where <?np will be referred to as the number of relativistic degrees of freedom (RDF). For BZ models we have 
gez ~ 100, both for the original example as well as others that have been studied using lattice Monte Carlo methods 0]. 
In the IR regime the RDF can be obtained using the AdS/CFT correspondence Q which gives 



gu = (7r 5 /8)(LM Pi ) 2 >100 



(7) 



since L, the AdS radius is > 2/Mp[. It follows that for all available models ^np^IOO, and in the IR we expect this 
number to be even larger. 



III. SM-NP INTERACTIONS; EQUILIBRIUM, FREEZE-OUT AND THAW-IN 

In order to understand the NP effects on cosmic evolution it is important to determine when and if this sector was 
in equilibrium with the SM sector. The standard approach for investigating this makes heavy use of the Boltzmann 
equation (BE) [6j, for our purposes the idea is to (i) calculate psm and pnp in terms of •& = Tnp — Ism using the BE 
and then (ii) use p oc T to obtain an evolution equation of the form 

t? + Am = -Ti9 (8) 

together with an explicit and calculable expression for the reaction rate T. H denotes the Hubble parameter, H 2 = 
[(87r/(3A/p ; )](psM + Pnp)j where we assumed a flat universe with zero cosmological constant. We will assume that 
the SM-NP interactions are of the form 



Ant — cCsmCnp 



(9) 
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The BE approach [6[ requires at intermediate steps the introduction of the unparticle distribution function, which 
might pose conceptual problems. This can be avoided by using instead the Kubo formalism Q, though it lacks the 
intuitive appeal of the BE. Both formalisms yield identical expressions for T (the momenta in the BE expression are 
defined in Fig. [5]): 



T — (— 1 \x- 

12T 4 VSSM flNP / 



iKubo = ds dt d d x ^O S M(-ia,x)OsM(t,0)^O NP (-i« J x)O N p(t,0: 

Zbe = \ d^ P d^ SM f3(E SM - E' SM ) 2 e-^ E \M\ 2 (2ir) 4 S(K - K') 



and one can show Tdu = T^..i... 
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(10) 



FIG. 2: Definition of momenta involved in the BE calculation of F. 

For the arguments below the detailed structure of T will not be needed and the following order of magnitude 
estimate suffices: 

p ^ e 2 A(ffSM + gNp) y 2 rfsM+2rf N p-7 ^ _ 9SM ffNP 
(^nsM+nNP-l <7SM5NP 

where ngM, «np ar e the number of SM and NP fields in C nlt , and 5g M , g^p are ^ ne RDF involved in the interaction 
(in the unparticle phase we take tinp = 2(4/ — 2, = du)- It then follows that 

r/# oc T 2dsM+2dNP ~ 9 (12) 

where the SM and NP sectors will be (de)coupled as long as T > H (T < H), the transition temperature Tf is 
determined by the condition T = H. Note that when g?sm + c?np < 4.5, the sectors are coupled for T below Tf, what 
we call a "thaw-in" scenario; in the complementary case dsM + <^np > 4.5 decoupling occurs for T < Tf (standard 
freeze-out scenario). 

For the calculation we will use 

( (6U)T„QQ/Mu BZ 

U ^ 



where Q denotes a quark in the BZ phase of the NP sector and On is the unparticle operator corresponding to 
QQ- Fot consistency we require Mu > Tf > An in the BZ sector, and An > Tf > v in the U sector (the lower 
value of v is needed because below this scale the SM operator <f>'<f> is no longer relevant). With these preliminaries 
one can determine the regions in the Mu — Ajj for which the SM and NP are coupled in each of the NP phases, the 
results are presented in Fig. [3] 



IV. UNPARTICLE EFFECTS IN BIG BANG NUCLEOSYNTHESIS (BBN) 

Let's first consider the case where SM and NP were in equilibrium down to a temperature Tf > v, and decoupled 
thereafter. The relationship between the NP and SM temperatures are thereafter determined by entropy conserva- 
tion [f|, specifically, 

(T f Rffg* NP (Tf) = (T NP R) 3 g* NP (T NP ) 
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FIG. 3: Regions in the Au — Mu plane corresponding to various freeze-out and thaw-in scenarios for du = 3/2, 2, 3, 7/2. Dark 
grey: SM-NP decoupling in the unparticle phase only; light gray: no SM-NP decoupling; in the white regions T u .f < v (Au, Mu 
are in TeV units). We assumed (?sm = 9bz = gu = 100, g$ M = 4, g' BZ — 50 and g' u = d U - BZ phase: n$Ai = n NP = 2, d SM = 2 
and d]vp = 3; hi phase: usm = 2, njvp = 2(du — 1), dsM = 2 and djvp = ciw- 



(T f R f y g * SM (Tf) = (T 7 i?) d 5 | M (T 7 ) 



(14) 



where gj^p and #g M stand for the NP and SM effective numbers of RDF conventionally p adopted for the entropy 
density Rf, (R) denote the scale factor at T = Tf (T 7 ), and T 7 , T^p denote the SM and NP temperatures at the 
BBN epoch. 

Using this we find 



T/vp — T 7 



37 9i + 9e + 9v 



1/3 



(15) 



_9i+9e 9sm{v) 

where g 7 = 2, g e = (7/8) x 4, = (7/8) x 3 x 2, (/sm(v) = 106.75. The NP contribution to the total energy density 
can be expressed in terms of an additional number of sterile neutrinos AN„ defined by the expression 
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Then, 



Current data 



9m 



9sm{v) 



4/3 



A TV,, 



.37 + 3e + 9u_ 

provides the limits AiV„ = ± 0.3 ± 0.3, leading to 

3iR < 20 



(16) 



(17) 



(18) 



at the 95% CL. In the extreme case where the SM and NP remain in equilibrium at the BBN epoch this bound is 
considerably strengthened: gm < 0.3. 

Viable unparticle models should exhibit conformal invariance with a small number of RDF in the IR. We are 
unaware of any model with these characteristics. In fact the AdS/CFT correspondence suggests this constraint is 
very strongly violated. 
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V. COMMENTS 



We have shown that strongly coupled NP can lead to T/H ~ T™ with n positive or negative, and that this results 
in a variety o freeze-out and thaw-in scenarios. Current BBN data generates strong constraints on the properties of 
the NP. Even for "normal" decoupling scenarios (n > 0) the BBN constraint is significant leading to a bound gm < 20 
to be compared to gm > 100 for the available models. 

Unparticle models also suffer from potential theoretical problems: the coupling to the SM necessarily breaks confor- 
mal invariance, but the scale at which this occurs lies below the BBN temperature for a range of du, so our conclusions 
still apply in this case. 

Acknowledgments 

We thank the Organizers of the XXXHI International Conference of Theoretical Physics for their worm hospitality 
during the meeting. This work is supported in part by the Ministry of Science and Higher Education (Poland) as re- 
search project N N202 006334 (2008-11) and by the U.S. Department of Energy grant No. DE-FG03-94ER40837. B.G. 
acknowledges support of the European Community within the Marie Curie Research & Training Networks: "HEP- 
TOOLS" (MRTN-CT-2006-035505), and " UniverseNet" (MRTN-CT-2006-035863), and through the Marie Curie Host 
Fellowships for the Transfer of Knowledge Project MTKD-CT-2005-029466. 



[1] H. Georgi, Phys. Rev. Lett. 98, 221601 (2007); Phys. Lett. B 650, 275 (2007). 

[2] T. Banks and A. Zaks, Nucl. Phys. B 196, 189 (1982). 

[3] J. C. Collins, A. Duncan and S. D. Joglekar, Phys. Rev. D 16, 438 (1977). 

[4] B. Svetitskv. larXiv:0901.2103l [hep-lat] . 

[5] S. S. Gubser, Phys. Rev. D 63, 084017 (2001) 

[6] E. W. Kolb and M. S. Turner, Addis on- Wesley (1990) 

[7] R. Kubo, J. Phys. Soc. Jap. 12, 570 (1957). R. Kubo, M. Yokota and S. Kakajima , J. Phys. Soc. Jap. 12, 1203 (1957). 

[8] G. Miele and O. Pisanti, Nucl. Phys. Proc. Suppl. 188 (2009) 15 larXiv:0811.4479l [astro-ph]]. 



